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The field of a eylindrieal dipole antenna in a compressible, isotropic, lossy plasma is represented 
by the field of a cylindrical current sheet immersed in the plasma. The antenna is short compared 
to a free space wavelength and for this reason a triangular current distribution is assumed. A formula 
for the input impedance is derived and compared with existing formulas for spherical, cylindrical 
and planar geometries. In addition, the effect of contact between the plasma and a metal antenna 
is estimated in the low frequency limit for the case of an antenna biased to the point of ion sheath 
collapse. 



1. Introduction 

Electroacoustic oscillations in isotropic, uniform plasmas arc characterized by the absence 
of an oscillating magnetic field. A propagating electroacoustic wave (plasma wave) is longitudinal 
and can exist at all frequencies above the plasma frequency if the electron gas is compressible 
(i.e., has a finite temperature). Electroacoustic effects have been studied extensively and a 
thorough discussion on the subject has been presented by Cohen [1961, 1962], whose work includes 
a calculation of the radiation resistance of a filamentary, sinusoidal current distribution; radiation 
resistance calculations have also been done by Chen [1963]. Hessel and Shmoys [1962] have 
calculated the field of an infinitesimal electric dipole in an infinite medium, and Hessel, Marcuvitz, 
and Shmoys [1962] have considered the problem of a magnetic current filament in free space over 
a compressible plasma half-space. Hall [1963] has obtained the impedance of a parallel plate 
capacitor and Fejer[1964] has carried out a similar computation for a single spherical electrode. 
Wait [1964] has studied the field of a slotted sphere and also [1965] the fields of finite and infinite 
cylindrical dipoles. In addition, Whale [1963], Mlodnosky and Garriott [1962], and Crawford and 
Mlodnosky [1964] have estimated the electroacoustic effect on impedance when an ion sheath is 
present. 

The small-signal impedance of an antenna in a plasma is determined principally by three 
phenomena, electroacoustic oscillations, electromagnetic oscillations and contact (or ion sheath) 
effects. At present it cannot be said that any one of these is negligible. In fact, all three pheno- 
mena would be closely coupled in any practical situation. Nevertheless, in the following analysis 
of a dipole antenna, contact phenomena will be disregarded at first and the antenna fields will be 
approximated by the fields of a given current distribution in a uniform medium. Under such 
conditions the electroacoustic and electromagnetic oscillations make separate and distinct contri- 
butions to the input impedance. 
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Research Laboratories under contract AF 19(628)-3900. Some of the material in this paper was contained in a dissertation submitted in partial fulfillment of the 
requirements for the Ph.D. degree, University of Illinois, Urbana, 111., 1963. 
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The approach outlined above requires a good estimate of the current distribution on the dipole. 
It is assumed that the dipole is short compared to a free space electromagnetic wavelength and that 
the electromagnetic oscillations alone determine the current distribution which, consequently, 
is approximately triangular (zero at the ends of the dipole and maximum at the center). Further- 
more, although the antennas considered here are short and thin compared to an electromagnetic 
wavelength, they are generally very long and moderately thick compared to an electroacoustic 
wavelength. Therefore, in order to take the thickness into account, the antenna must be repre- 
sented by a cylindrical sheet of current and not by a current filament. 

2. Theoretical Development 

In rationalized MKS units and for e> w time variation, Maxwell's equations are 

Vx£=- JCD/MoH (1) 

\7xH=j<x)eoE-Nev + J (2) 

in which J is the source current, TTthe electron velocity (the ions are stationary), e the electron 
charge magnitude and TV the average electron density in the plasma. Conservation of the number 
of electrons (continuity) is expressed by 

NV •v+ja>n = (3) 

in which n is the oscillating part of the electron density. The equation of motion is 

ja)Nmv = — NeE — Vp — Nmvv (4) 

in which p is the scalar pressure, m the electron mass and v the collision frequency. The adia- 
batic equation of state is 

p = ykTn (5) 

where k is Boltzmann's constant and T is the electron temperature. The factor y is the ratio of 
specific heats and may be set equal to 3 for one-dimensional adiabatic compression [Spitzer, 
1962]. Combining the equations of motion and state gives 

jo)UNmv = - NeE- mV^Vn (6) 

in which U=l—jZ, Z = — , and V 2 = . The quantities E and v may be eliminated from (2), 

a) m 

(3), and (6) to yield the differential equation for n, 

(V 2 -a 2 )n = -^r i V-J 

jcoeomV 2 (7) 

in which a* = — ,X=Jsmdo> i N = m - 

At this point it is convenient to separate E and v into electromagnetic and plasma (electro- 
acoustic) components, designated by the subscripts e and p, respectively [Cohen, 1961]. Thus (6) 
can be written as two equations: 

jcoUNmvo = — NeE e (8) 

jcoUNmVp = — NeE p — mf^Vn . (9) 

560 



The differential equations now may be separated into two groups. The electromagnetic group is 

VXEe=-j(0/JLoH (10) 

VXfl =j<D€<>KoEe + J (11) 

and the electroacoustic group consists of (7) and 

V 2 

^UK^N Vn (12) 



v P =J 



- mV 2 (l-K ) 

Ep= NeK V " 



(13) 



in which K =l— XU ! . Combining (7) and (13) gives the differential equation for the electro- 
acoustic part of the electric field strength: 



7&>e A 



(14) 



In order to evaluate the electroacoustic contribution to the impedance of a dipole it is neces- 
sary to select an appropriate current density function J, find the component of E p parallel to J 
by solving (14) and evaluate the impedance using the Poynting theorem. If J is taken to be in 
the z-direction and if E z is the z component of E p , then (14) becomes 



(V 2 -a 2 )E z 



1-XodV* 



jco€oK dz 2 
For the case of unit input current, the current density is given by 



(15) 



'-^'M 



(16) 



where p is the dipole radius, r is the cylindrical radial coordinate and J(z) is the triangular func- 
tion shown in figure 1. 

Equation (15) will be solved using the transform pair 



f(7,k) = j X j~f(r,z)e-^J (yr)rdrdz 



f(y, k)e> k:: Jo(yr)ydydk. 



(17) 
(18) 




FIGURE 1. The assumed longitudinal current distribution 
function. 
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The transformed current density is 

^ = 2^zfe (e ~ JAL4 "^ L " 2) ' /o(7p) < 19 ) 

in which Jo(yp) is the Bessel function of the first kind. The electric field can be expressed as 



K —l ( x ( x e~ jkL -h e jkL — 2 

E * = ~ v t<o \2 ^2-l 2-u 2 & kz Myp)Jo(yr) ydydk. 

jaj€oKoL(Z7rr J -x Jo k* + y l + a l 



(20) 



In order to simplify the calculations, the impedance of a monopole of height L over a ground plane 
will be computed and the dipole impedance will be obtained by doubling the monopole impedance. 
Thus the field calculations may be limited to the range ^ z ^ L. Furthermore, for impedance 
calculations only the field at r = p is needed. If, in addition, Vy 2 + a 2 is taken to have a positive 
real part then integration of (20) with respect to k gives 



/£ ft — 1 f «= e -(L-z)\ ^2+ a 2 _[_ e _(L+2)Vy2 +a 2 _ 2 e -*Vy2+a2 



The use of the formula 



2 r 

Jl(yp)=-\ M2ypcos0)d0 (22) 

rr Jo 

and the assumption that a has a positive real part permits integration with respect to y (it is a form 
of "Sommerfeld's integral"): 



n - 2 p -a \/(2r-L)2+(2p cos 0)2 

/(L)=2f -d». 



^here 



ttJo V(z-L) 2 +(2p cos 0) 2 

Integration with respect to 6 will be delayed in order to simplify subsequent calculations. 

For the case of unit input current the dipole impedance due to electroacoustic oscillations 
is found by integrating the product {—J*E Z ) over the volume occupied by the current, a procedure 
which may be deduced from the Poynting theorem. Thus 



Z?=~2 



j/ l-j\E z dz. (24) 



L 



Under the assumption that L 2 > p 2 , integration with respect to z gives 



Zp = 



Ko-l 2 



j<t)7T€oK { 



[_2 r 

oLtt J 



K (2ap cos 0) - 2E x {aL) + E&aL) 

2aL 



. Q e -2apcos$ 



d0. (25) 



In this formula the modified Bessel function of the second kind is given by 

Too e~ au 

K ^ x n x v7^ du (26) 
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and the exponential integral is given by 



Et(aL) 



J L U 



du. 



(27) 



Integration with respeet to 6 gives 
Ko-1 



Zt>=: 



jcoireoKoL 



Io(ap)K (ap) - 2Ei(aL) + E,(2aL) 



+ 2~L 4e ~ aL " e ~ 2aL ~ 3I °( 2a P) + 3L (2ap) 



(28) 



in which I {x) is a modified Bessel function of the first kind and Lq(x) is a modified Struve function. 
In most cases of interest ah is very large and ap is fairly small; under such conditions only the 
first term in (28) is significant. 

A complete impedance formula can be obtained by combining the first term of (28) with the 
electromagnetic impedance contribution Z m . For the case of a very short dipole, Zf n is pre- 
dominantly reactive and the expression for it may be deduced readily from the free space expression 
[see for instance Schelkunoff and Friis, 1952]. Thus the complete formula is 



7 = 7 e 4- 7 P ■ 



I 



J(07T€qKoL 



In \-(\-KMap)K<Aap) 

P 



(29) 



This expression is convenient to use when o> < a>#. However, when a> > a>\, it is preferable to set 
a=jp. Thus (29) becomes 



7j\n • 



1 



jojTreoKoL 



/7l--l+(l 

P 



-Xo)|{yo08p)yVo08p)+7yS08p)} 



(30) 



For a lossless plasma, Z xn has a positive real part which may be associated with the radiation of 
electroacoustic waves. 

Similar formulas have been worked out for a lossless medium by Hall [1963] and Fejer [1964] 
who used boundary-value techniques to derive impedance formulas for a single sphere and a par- 
allel-plate capacitor, respectively. The following is a summary of the various formulas; note 
that R is sphere radius and S and D are parallel plate area and spacing, respectively. 

X 



Sphere: Zi n = -r 
Plates: Z\ n 



_1 

jO)^7T€o(l~X)R 

2D 



Dipole: Zj n = ~ 



jo)6 (l—X)S 

1 



1 



1+aR 

X tanh aD 
Trf) 



JO)TT€o{ 1 — X)L 



In 1 — XIo(ap)Ko(ap) 



(31) 
(32) 

(33) 



Of the three expressions only (32) has no positive real part for a> > a) N . The similarity of the 
formulas is most apparent if to < a>\ and if the temperature approaches zero {a—* °°): 



Sphere: Z\ n =-. 
Plates: Z xn 



1 



j(0^7T€o(l-X)R 

2D 



aR 



X 



Dipole: Z m = 



ja)€o(l— X)S\_ aD 
1 
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ln--\ 



X 

2ap 



(34) 
(35) 
(36) 




2.5 3.0 




FIGURE 2. Impedance calculations for N = 10 4 el. /cm. 3 



FIGURE 3. Impedance calculations for N = 10 6 eljcm. 3 



3. Discussion of Results 

Equation (36) is helpful in estimating the magnitude of electroacoustic effects; for instance, 
if X is set equal to 2, a value of {ocp)~ l of the order of unity clearly indicates an appreciable effect. 



{ap) w = 2^p 



at X = 2. 



(37) 



Representative parameters in the ionosphere are T= 300 °K, X = 300 m and p = 1 cm. Under these 
conditions, (ap)~ 1 = 1.86 and thus strong electroacoustic effects may be anticipated. In labora- 
tory plasmas [see, for instance, Balmain, 1964] representative parameters are r=300°K, \ = 30 
cm, and p = l mm. Under these conditions, (ap) _1 = 0.0186 and thus electroacoustic effects are 
negligible according to the theory: developed so far. 

Computations were carried out for a dipole with dimensions p = 5 mm. and I = 5m. and the 
results are shown in figures 2 and 3. In the computations ^ = and TV and T are representative 
of the ionosphere. Figure 2 shows a low frequency resonance at X ~ 2 (all the curves would indi- 
cate a resonance if the X axis were extended far enough). Such a resonance has been noticed 
by Fejer [1963] and is contained in the formula derived by Hall [1963]. Inspection of (31) to (36) 
reveals that the resonant frequency depends on the size and shape of the antenna. 

The radiation of electromagnetic waves has not been included in the discussion so far but its 
effect can be estimated by the use of the familiar formula for the radiation resistance of a short 
dipole [Jordan, 1950, for instance]. Upon insertion of the relative permittivity of a lossless plasma, 
the radiation resistance formula becomes 



D 807T 2 L 2 Vr=X 

R — \> 



(38) 



in which X is the free space wavelength. With suitable normalization (38) is plotted on figure 3 
where it can be compared directly with the function G. The effect of electromagnetic radiation 
is appreciable only for small X in figure 3 and is entirely negligible in figure 2. 

It is worthwhile to compare the electroacoustic radiation resistance as given by (30) with the 
resistances calculated by other authors for filamentary source currents. Cohen [1962], Chen [1963] 
and Wait [1965] all have carried out such calculations, and in all three cases the electroacoustic 
radiation resistance approaches a finite value as the electron temperature approaches zero. This 
behavior is entirely the result of having considered a source current of infinitesimal radius rather 
than one of finite radius. Equation (30) indicates that, for a dipole of finite radius, the entire electro- 
acoustic part of the input impedance approaches zero as the electron temperature approaches 
zero. For finite temperature and vanishingly small radius, however, the electroacoustic radiation 
resistance in (30) is identical to that derived by Wait [1965] for a filamentary source current. 
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The computations shown in figures 2 and 3 indicate that the input reactance is very strongly 
affected by the compressibility of the plasma. This theoretical result is ample cause for a re- 
examination of the assumptions made at the beginning of the analysis. The assumption of a 
triangular current distribution is particularly suspect since it implies that the electroacoustic 
oscillations do not affect the current distribution. Without electroacoustic effects, the triangular 
current assumption is justified provided that the antenna is short compared to a wavelength under 
propagating conditions and short compared to the penetration length under cutoff conditions. 
However, the strong electroacoustic effect in the impedance suggests that there may be an equally 
strong electroacoustic effect in the current distribution. Such a conclusion has also been reached 
by Wait [1965] following his analysis of an infinitely long dipole. 

Also open to question is the assumption that the antenna can be represented by a current 
sheet in a uniform medium. This is the approach used in the "induced EMF" method of impedance 
calculation, and such an approach to a certain extent neglects the physical presence of the antenna 
surface; consequently a high order of accuracy cannot be expected in the reactance calculations. 
Furthermore, the equation of continuity shows that the antenna model is one in which the electrons 
in the medium are not taken into the metal surface to become part of the antenna current. In 
this respect the model is similar to the one used by many authors who employ the rigid boundary 
condition [Fejer, 1964, and Hall, 1963 for instance]. It is evident that a better plasma-metal 
boundary condition must be found before full confidence can be placed in any impedance calculation. 

4. Contact Effects 

Any antenna immersed in a plasma is surrounded by an ion sheath and consequently the 
electron density is highly nonuniform within a few Debye lengths from the surface. If a uniform- 
medium impedance theory is to be compared with experiment, the sheath must be collapsed by 
giving the probe a positive d-c bias with respect to some large reference electrode. Although 
sheath collapse ideally results in a uniform electron density adjacent to the antenna surface, 
it also introduces an appreciable steady flow of electron current. If the plasma is compressible 
the application of an RF signal causes a density modulation and consequently an RF current in 
the stream of electrons flowing to the antenna surface. 

The impedance contribution due to this effect may be computed very easily in the low fre- 
quency limit. When (o = Q and v = Q, (6) becomes 

— NeE = mV 2 Vn. (39) 

If E =— V where i// is a scalar potential, integration of (39) between two plates of a capacitor gives 

/Ve(i//i — 1// 2 ) = mV^irti — m). (40) 

For a symmetrical geometry it is possible to set 1//2 — — 1//1 and ni — —n\ from which 

Neifji = mV 2 ni. (41) 

The oscillating current to an electrode with area A may be expressed as 



VkT 
2™ («) 



Thus the input resistance for a cylindrical dipole is given by 

_2\lj 1 _ 2yV2TrmkT _ 2y fmkT 
Kin ~ I " ANe 2 ~pLNe 2 V 2tt ' < 43 ) 



/ ANe 2 pLNe 2 V 2tt 
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Rin would constitute the total input impedance in the limit as frequency approaches zero and at 
higher frequencies one would expect to find a shunt resistance similar to Ri n . The influence of the 
high frequency resistance is difficult to predict but, since Ri n a /V _1 , little contact effect would be 
expected at low electron densities (say 10 4 el/cm 3 ). 

Under isothermal conditions y is equal to unity and thus (43) reduces to the resistance dis- 
cussed by Mlodnosky and Garriott [1962 J whose formula is derived by taking the slope of the Lang- 
muir probe voltage-current characteristic. Thus one might expect a transition from adiabatic to 
isothermal behavior as frequency is lowered; calculations by Pavkovich [1964] suggest that there 
is such a transition but that it is very gradual as frequency changes from (On to zero. It should be 
noted that the isothermal resistance has been used by Balmain [1964] to account for high apparent 
losses in laboratory impedance probe experiments. Evidently, contact effects are appreciable 
even in laboratory plasmas, and furthermore these effects are closely related to electroacoustic 
oscillations. 
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